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Chandrasekhar separated the Dirac equation for spinning and massive particles in Kerr geometry 
in radial and angular parts. Chakrabarti solved the angular equation and found the correspond- 
ing eigenvalues for different Kerr parameters. The radial equations were solved asymptotically by 
Chandrasekhar. In the present paper, we use the WKB approximation to solve the spatially com- 
plete radial equation and calculate analytical expressions of radial wave functions for a set of Kerr 
and wave parameters. From these solutions we obtain local values of reflection and transmission 
coefficients. 
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I. INTRODUCTION 



One of the most important solutions of Einstein's equation is that of the spacetime around and inside an isolated 
black hole. The spacetime at a large distance is flat and Minkowskian where usual quantum mechanics is applicable, 
while the spacetime closer to the singularity is so curved that no satisfactory quantum field theory could be developed 
as yet. An intermediate situation arises when a weak perturbation (due to, say, gravitational, electromagnetic or 
Dirac waves) originating from infinity impinges on a black hole, interacting with it. The resulting wave is partially 
transmitted into the black hole through the horizon and partially scatters off to infinity. In the linearized ('test field') 
approximation this problem has been attacked in the past by several authors [1-4]. The master equations of Tcukolsky 
[2] which govern these linear perturbations for integral spin (e.g., gravitational and electromagnetic) fields were solved 
numerically by Press & Teukolsky [5] and Teukolsky & Press [6]. While the equations governing the massive Dirac 
particles were separated by Chandrasekhar [3]. So far, only the angular eigenfunction and eigenvalue (which happens 
to be the separation constant) have been obtained [7]. Particularly interesting is the fact that whereas gravitational 
and electromagnetic radiations were found to be amplified in some range of incoming frequencies, Chandrasekhar [4] 
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predicted that no such amplifications should take place for Dirac waves because of the very nature of the potential 
experienced by the incoming fields. However, these later conclusions were drawn using asymptotic solutions and no 
attempt has so far been made to determine the nature of the radial wave functions, both incoming and outgoing, for 
the Dirac wave perturbations. He also speculated that one needs to look into the problem for negative eigenvalues 
(~X) where one might come across super-radiance for Dirac waves. 

In the present paper, we revisit this important problem to study the nature of the radial wave functions as a function 
of the Kerr parameter, rest mass and frequency of incoming particle. We also verify that super-radiance is indeed 
absent for the Dirac field. Unlike the works of Press & Tcukolsky [5] and Teukolsky & Press [6] where numerical 
(shooting) methods were used to solve the master equations governing gravitational and electromagnetic waves, we 
use an approximate analytical method for the massive Dirac wave. The details of the method would be presented 
below. 

The plan of the paper is as follows: in the next Section, we present the equation governing the Dirac waves (waves 
for half-integral massive spin particles) as they were separated into radial and angular co-ordinates. We then briefly 
present the nature of the angular eigenvalues and eigenfunctions. In §3, we present our method of solution and present 
the spatially complete radial wave functions. Finally, in §4, we draw our conclusions. 



Chandrasekhar [3] separated the Dirac equation in Kerr geometry into radial (R) and angular (S) wave functions. 
Below, we present these equations from Chandrasekhar [4] using the same choice of units: we choose H = 1 = G = c. 
The equations governing the radial wave- functions R±i corresponding to spin ±| respectively are given by: 



II. THE DIRAC EQUATION IN KERR GEOMETRY 



A5p i?_i = {X + im p r)A^R + i, 



(la) 



A^vlA^R + i = (X - im p r)R_i, 



(lb) 



where, the operators V n and T>\ are given by, 



T> n = d r + — + 2n 



(r - M) 
A 



(2a) 
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and 

A = r 2 + a 2 - 2Mr, (3a) 

K = (r 2 + a 2 )a + am. (36) 

Here, a is the Kerr parameter, n is an integer or half integer, a is the frequency of incident wave, M is the mass of the 
black hole, m p is the rest mass of the Dirac particle, A is the eigenvalue of the Dirac equation and m is the azimuthal 
quantum number. 

The equations governing the angular wave- functions S± 1 corresponding to spin ±| respectively are given by: 

ClS + i = — (A — am p cos 6)S_ i (4a) 

£t S_i = +(A + am„ cos 0)S+ i (46) 

2 2 '2 

where, the operators C n and £jj are given by, 

C n = do + Q + n cot 0, (5a) 

Ct = de-Q + ncot9 (56) 

and 

Q = aa sin + m cosec 6. (6) 

Note that both the radial and the angular sets of equations i.e., eqs. l(a-b) and eqs. 4(a-b) are coupled equations. 
Combining eqs. 4(a-b), one obtains the angular eigenvalue equations for the spin-i particles as [7] 



t am p sm8 f 2 2 2 2 

£±£i + -C\ + (A - a z mt, cos z 6) 

2 2 A + am„ cos & 2 F 



5_i=0. (7) 



There are exact solutions of this equation for the eigenvalues ^ and the eigenfunctions S_ 1 when p = ^ = 1 in 
terms of the orbital quantum number I and azimuthal quantum number m. These solutions arc [7]: 

1, 



i 



2(1 + 1) + aax 



(8) 



and 



l5lm = iYl m - — — — — iVj+lm (9) 

2 2 2(1 + 1) + aax 2 

where, 

p = F(l,l); x = F(l + 1,1 + 1); y = F(l 1 l + l) 

and 

F(/!,Z 2 ) = \(2l 2 + l)(2h + 1)]* < Z 2 l™0|?im > 

[< ^1^0|ii^ > +(-l) i2 - ; < i 2 lmO|iim > [< hl\®\h\ > +(-l) h - l P V2 < l 2 l - >]]. (10) 

with < ....|.. > are the usual Clebsh-Gordon coefficients. For other values of p one has to use perturbation theories. 
Solutions upto sixth order using perturbation parameter aa is given in Chakrabarti [7]. The eigenfunctions X are 
required to solve the radial equations which we do now. 

This radial equations l(a-b) are in coupled form. One can decouple them and express the equation either in terms 
of spin up or spin down wave functions R±l but the expression loses its transparency. It is thus advisable to use the 
approach of Chandrasekhar [4] by changing the basis and independent variable r to 

2Mr + + am/a, (r \ 2Mr-+am/a, (r \ . . 

r* = r H log 1 log 1 , (11 

r + — r- \r+ J r + — r- \r_ / 



(for r > r+), 



_d_ _ A__d 

dr* lo 1 dr ' 



(12) 



LU 2 = r 2 + a 2 ; a 2 = a 2 + am/cr, (13) 
to transform the set of coupled equations 1 (a-b) into two independent one dimensional wave equations given by: 



-f— + ia] P_i = ^-(X + im p r)P,i. 
dr* I 2 w 2 v p i + 2 



(15) 



Here, 2?o = ^(377 + * cr ) an( l ^0 = 7r(d?7 — * cr ) were use d an d wave functions were redefined as R_i = P_i and 

Asp, 1 = P, 1. 

^2 ^2 

We now define a new variable, 



tan 1 (m p r/X) 



(16) 



which yields, 



costf = 



and sin 9 = 



and 



(A~ ± zm p r) = exp(±z6')v / (^ + m_r ), 



(17) 



so the coupled equations take the form, 



- P +i = ^J(A 2 + my) 1/2 ^i«P 
ax* / +2 w 2 p 2 



-i tan 



_i / '"'pi 



(18a) 



and 



i<A P_i = (A 2 + m 2 r 2 ) 1 / 2 P + i exp 



rfr 



i tan 



1 / m p r 



(186) 



Then defining, 



P + i = V+i exp 



1 . , _i / w„r 
--z tan 



(19a) 



and 



P_i = V'-i exp 



1 1 / w„r 

+-, tan- f ^ 



(196) 



we obtain, 



dr* \ 



2cr A 2 + m 2 r 2 



(20a) 



and 



w 2 2ct A 2 + m 2 r 2 



A* 



(206) 
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Further choosing f* = r* + 7^-tan 1 (Hh^L\ so that dr* = (1 + _a i )rfr«, the above equations become, 



■^r-iv)il> + i=W1>_i, (21a) 



and 



where, 



dh +ia)i,_ h =W^ +h . (216) 



Al(^ + m ^)3/ 2 
w 2 (* 2 + m 2 r 2 )+*TO p A/2a V ' 



Now letting Z± = ip + i ± we can combine the differential equations to give, 



W) Z+ = iaZ-, (23a) 



df 



and 



4- + W) Z_=i(jZ + . (236) 
^From these equations, we readily obtain a pair of independent one-dimensional wave equations, 



where, V± = W 2 ± dw 



A^fi 2 + m 2 r 2 ) 3 / 2 
[w 2 (* 2 +m 2 r 2 ) +}to p A/2ct] 



— - + a 2 )Z± = ViZ±. (24) 

GET** 



[A^(} 2 + m 2 r 2 ) 3 / 2 ± ((r - M){% 2 + m 2 r 2 ) + 3m 2 rA)] 



A§0\ 2 + m 2 r 2 ) 5/2 

V^+W^W ^ + ^ + ^ + ^ - M) M (25) 
One important point to note: the transformation of spatial co-ordinate r to r* (and f*) is taken not only for 
mathematical simplicity but also for a physical significance. When r is chosen as the radial co-ordinate, the decoupled 
equations for independent waves show diverging behaviour. However, by transforming those in terms of (and f*) 



we obtain well behaved functions. The horizon is shifted from r = r + to f* = — oo unless a < a s = —am/2Mr + (eq. 
11). In this connection, it is customary to define a c where a 2 — (eq. 13). Thus, a c = —raja. If a < a s , the region 
is expected to be super-radiant [4] because for integral spin particles for a < sigma s there exhibit super-radiation. 

III. SOLUTION OF THE RADIAL EQUATION 

Out of the total physical parameter space, in one region (region I) the total energy of the particle is always greater 
than the height of potential barrier and in the other region (region II) the energy is less than of the maximum height 
of the potential barrier. In region II, the wave hits the wall of barrier and tunnels through it. One has to treat these 
two cases a little differently. 

The usual WKB approximation [8] is used to obtain the zeroth order solution. We improve the solution by properly 
incorporating the inner and outer boundary conditions. After establishing the general solution, we present here the 
solution of eq. (4) for three sets of parameters as illustrative examples. For those examples the choice of parameters 
is made in such a way that there is a significant interaction between the particle and the black hole, i.e., when the 
Compton wavelength of the incoming wave is of the same order as the radius of the outer horizon of the Kerr black 
hole. So, 

G[M + 7(M 2 - a 2 )] n 



rripC 



(26) 



We choose as before G = H = c = 1, so 



m p F . (27) 

P [M+^[M 2 -a 2 )] 

Similarly, the frequency of the incoming particle (or wave) should be of the same order as the inverse of the light 
crossing time of the radius of the black hole, i.e., 

„3 

a. (28) 



G[M + y/(M 2 - a 2 ) 
Using the same units as before, we can write, 



m 



v ~ o ~ [M + V(M 2 - a 2 )}- 1 . (29) 



In principle, however, one can choose any values of a and m p for a particular black hole and the corresponding solution 
is possible. 



One can easily check from equation (25) that for r — > oo (i.e., f* — > oo) V± — ► m^. So we expand the total parameter 
space in terms of the frequency of the particle (or wave), a and the rest mass of the particle, m p . It is clear that in 
half of the parameter space spanned by a — m p where, a < m p , particles are released at finite distance with so little 
energy that they cannot escape to infinity. In this case, the total energy <~ er 2 of the incoming particle at a large 
distance is less than the potential energy of the system. We will not discuss solutions in this region. The rest of the 
parameter space [a > m p ) is divided into two regions - I: E > V m and II: E < V m , where E is the total energy of 
the incoming particle and V m is the maximum of the potential. In Region I, the wave is locally sinusoidal because the 
wave number k is real for the entire range of f*. In Region II, on the other hand, the wave is decaying in some region 
when E < V, i.e., where the wave 'hits' the potential barrier and in the rest of the region, the wave is propagating. 
We shall show solutions in these two regions separately. In Region-I whatever be the physical parameters, the energy 
of the particle is always greater than the potential energy and the WKB approximation is generally valid in the whole 
range (i.e. p- Jr- << 1)- In cases of Region-II, the energy of the particle is always less than the maximum height 
of potential barrier. Thus, at two points (where, k = 0) the total energy matches the potential energy and in the 
neighbourhood of those two points the WKB approximate method is not valid. They have to be dealt with separately. 
In Fig. 1, we show contours of constant w max =max(p^) for a given set (<r, m p ) of parameters. The labels show 
the actual values of w max . Clearly, in most of the parameter regions the WKB approximation is safely valid for any 
value of tv One has to employ a different method (such as using Airy Functions, see below) to find solutions in those 
regions where w ma x attains a large value which indicates the non-validity of WKB method. 

Solutions of Region I 

We re-write equation (24) as, 

^± + (a 2 - V + )Z + = 0. (30) 

This is nothing but the Schrodinger equation with total energy of the wave ct 2 . This can be solved by regular WKB 
method. 

Let fc(r») = y/(a 2 — V+), u(r») = J fc(f*)df* + constant, k is the wavenumber of the incoming wave and u as the 
Eikonal. The solution of the equation (30) is, 

Z + = ^-exp(iu) H —exp(-iu). (31) 
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with 

A\ + A 2 _ = k. (32) 

The motivation of equation (32) is to impose the WKB method at the each space point so that sum of the transmission 
and reflection coefficients are same at each location. In this case a 2 > V + all along and also \ « k, so the WKB 
approximation is valid in the whole region. 

It is clear that a standard WKB solution where A + and A- are kept constant throughout, can not be accurate in 
whole range of r» , since the physical inner boundary condition on the horizon must be that the reflected component is 
negligible there (since there the potential barrier height goes down to zero). Thus the WKB approximation requires 
a slight modification in which a spatial dependence of A± is allowed. On the other hand, at a large distance, where 
the WKB is strictly valid, A + and A_ should tend to be constants, and hence their difference is also a constant: 

A+ - A- = c. (33) 

Here, one can choose also the sum of A + and A_ are constant instead of difference as equation (33), but the final 
result will not be affected. Here, c is determined from the WKB solution at a large distance. For simplicity we choose 
A±s are real. This along with (32) gives, 

4(r) = ± £ + « (34) 

This spatial variation, strictly valid at large distances only, should not be extendible to the horizon without correcting 
for the inner boundary condition. These values are to be shifted by, say, A±h respectively, so that on the horizon one 
obtains the physical R and T. We first correct the reflection coefficient on the horizon as follows: Let A_h be the 
value of A_ on the horizon (see, eq. 34), 



c y/[2fc(r + ) - t?] 
A - h = -2 + 2 ' 

It is appropriate to use A- = A_ — A_ h rather than A_ since A- vanishes at r = r+. 

Incorporating these conditions, the solution (31) becomes, 

A+ A- 

Z + = —-exp(iu) H —exp(-iu). (35) 

with the usual normalization condition 

A\ + A 2 _= g. (36) 



where, A+ = A+ - A +h . 

Determination of A+h is done by enforcing R obtained from eq. (37a), which is shown below, is the same as that 
obtained by the actual WKB method. The q is used to compute the transmission coefficient T from eq. (36). In this 
way, normalization of R + T = 1 is assured. 

The normalization factor q — > k as f* — > oo and the condition << q is found to be satisfied whenever 

\w~ « k 1S satisfied. This is the essence of our modification of the WKB. In a true WKB, A± are constants and the 
normalization is with respect to a (almost) constant k. However, we are using it as if the WKB is instantaneously valid 
everywhere. Our method may therefore be called 'Instantaneous' WKB approximation or IWKB for short. Using the 
new notations, the instantaneous values (i.e., local values) of the reflection and transmission coefficients are given by 
(see, eq. 35), 



The variation of reflection and transmission coefficients would be well understood if we imagine the potential barrier 
consists of a large number of steps. From simple quantum mechanics, in between each two steps, we can calculate the 
reflection and transmission coefficients [9]. Clearly these reflection and transmission coefficients at different junctions 
will be different. This is discussed in detail below. To be concrete, we choose one set of parameters from Region I. 
Here, the total energy of the incoming particle is greater than the potential barrier height for all values of f*. We use, 
Kerr parameter, a — 0.5; mass of the black hole, M = 1; Mass of the particle, m p = 0.8; orbital angular momentum 
quantum number, I = |; azimuthal quantum number, m = — \\ frequency of the incoming wave, a — 0.8. The derived 
parameters are, r+ = M + ^(M 2 - a 2 ) = 1.86603; cr c = 1; a s = 0.066987; a 2 = -0.0625. For these parameters, the 
eigenvalue is ^ = 0.92 [7]. 

Here it is clear that a is in between a c and a s and a 2 < 0, r+ > \a\. So we are in a strictly non super-radiant 
regime since here, a > <r s [4] . 

,/From eq. (24) we observe that there are two wave equations for two potentials V+ and V-. The nature of the 
potentials is shown in Fig. 2a. It is clear from the Fig. 2a that the potentials V± are well behaved. They are 
monotonically decreasing as the particle approaches the black hole, and the total energy chosen in this case (a 2 ) is 




(37a) 




(376) 



Whatever may be the value of the physical parameters, \-%jh « k is satisfied in whole range of f* for region I. 
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always higher than V±. For concreteness, we study the equation with potential V+. A similar procedure (IWKB 
method) as explained above can be adopted using the potential V- to compute Z_ and its form would be 

A' — A' A' — A' 

Z_ = + +h exp(iu') ~ ~ h exp{-iu'). (35') 

vQ vl 

Note the occur ance of the negative sign in front of the reflected wave. This is to satisfy the asymptotic property of 
the wave functions. 

In Fig. 2b, we show the nature of V + (solid curve), k (dashed curve) and E{= a 2 ) (short-dashed curve). In the 
solutions (eq. 35 and 35') the first term corresponds to the incident wave and the second term corresponds to the 
reflected wave. 

In Fig. 2c, the variation of reflection and transmission coefficients are shown. It is seen that as matter comes close 
to the black hole, the barrier height goes down. As a result, the penetration probability increases, causing the rise of 
the transmission coefficient. 

Local values of the reflection and transmission co-efficients could also be calculated using the well known quantum 
mechanical approach. First one has to replace the potentials (as shown in Fig. 2a) by a collection of step functions 
as shown in Fig. 3a. The standard junction conditions of the type, 

Z+, n = Z +>n+ \ (38a) 



where, 



^+,n — A n exp[ik n r* n ] -t- J3 n exp[ ik n v*_ n \ 



and 



dZ+ _dZ+ 

dfH"-^H" +1 {386) 



where, 



dZ 

77 \n — ik n A n G.xp(ik n T* : n ) ik n B n €xp( ik n v* n ) 

at each of the n steps were used to connect solutions at successive steps. ^From the simple quantum mechanical 
calculation we obtain the reflection and transmission coefficients at the each junctions. Clearly at different junctions 
i.e., at different radii this reflection and transmission coefficients will be different. As before, we use the inner boundary 

II 



condition, to be R — > at f* — > — oo. In reality we use as many steps as possible to follow accurately the shape of the 
potential. Smaller step sizes were used whenever k varies faster. Fig. 3b shows the comparison of the instantaneous 
reflection coefficients in both the methods. The agreement shows that the WKB can be used at each point quite 
successfully. 

It is to be noted, that, strictly speaking, the terms 'reflection' and 'transmission' coefficients are traditionally 
defined with respect to the asymptotic values. The spatial dependence that we show are just the dependence of the 
instantaneous values. This is consistent with the spirit of IWKB approximation that we are using. 

The radial wave functions R + ± and R_i which are of spin up and spin down particles respectively of the original 

Dirac equation are given below, 

i a + cos(u - 9) + a_cos(w + 9) a' + cos(V - 9) - a'_cos(u' + 9) 
Re(R^) = ^= + ^= (39a) 

i a + sin(u - 9) - a_sin(w + 9) a' + sin(V - 9) + a'_sm(u' + 9) 
Im(R h A>) = ^ + -j= (396) 

a + cos ( u + 0) + a - cos ( u - 9) a+cos(w' +9) — a'_cos(u' - 9) 
MR -i ] = U% 27P (39C) 

a + sin(w + 9) - a_sin(w - 9) a+sin(u' + 9) + a'_sm(u' - 9) 

/m(i? - ) = 271 VF m) 

Here, a+ = (A + - A +h )/ \/{q/k) and a_ = (A_ - A- h )/ ^f{q/k). -^t and ^= are the transmitted and reflected 
amplitudes respectively for the wave of corresponding potential V- . 

In Fig. 4(a-d) we show the nature of these wavefunctions. The eikonals used in plotting these functions (see, 
eq. 39[a-d]) have been calculated by approximating V± in terms of a polynomial and using the definition u(f*) = 
fy/l<r*-V±)dr*. This was done since V± is not directly intcgrable. Note that the amplitude as well as wavelength 
remain constants in regions where k is also constant. 

Solutions of Region II 

Here we study the solution of a region where for any set of physical parameters, the total energy of the incoming 
particle is less than the maximum height of the potential barrier. So the WKB approximation (more precisely, our 
IWKB approximation) is not valid in the whole range of f* . In the regions where the WKB is not valid, the solutions 
will be the linear combination of Airy functions because the potential is a linear function of f* in those intervals. At 
the junctions one has to match the solutions including Airy functions with the solution obtained by WKB method. 

12 



In the region where the WKB approximation is valid, local values of reflection and transmission coefficients and 
the wave function can be calculated easily by following the same method described in previous sub-section (solutions 
of region I) and the solution will be same as equation (35, 35'). In other regions, the equation reduces to 

<,2Z " xZ + = (40) 



df 2 

where x = f3^(r*—p), (3 is chosen to be positive and p is the critical point where the total energy and potential energy 
are matching. 

Let Z + (x) = x?Y(x) and considering region x > the equation (40) reduces to 

* 2 S + ^-f* 3 + 7V(-) = 0. (41) 



dx 2 dx \ 4 

By making another transformation 



i=\xi (42) 



we obtain 



d 2 Y AY /, 1 



? w +^-{e+s)nz) = o, (43) 

this is the modified Bessel equation. The solution of this equation is and So the solution of cq. (40) 

will be 

Z + (x)=x*[C 1 I + i(0 + C 2 I_i{Q]. (44) 
When x < the corresponding equation is, 

?w +(d i + { (2 -lh () =°- (45 » 

which is the Bessel equation. The corresponding solution is 

Z + (x) = M*[AJ + i(0+£>2./_i(0], (46) 

where J ± i and J±g are the Bessel functions and the modified Bessel functions of order | respectively. 
The Airy functions are defined as 

Ai{x) = ±xt[I_i{0-I + i(0], x>0, (47) 
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Ai{x) = ±\x\t[J_i(t) + J + i(t)], x<0, (48) 



Bi(x) = ±x±[I_i(t) + I + i(t)], x>0, (49) 

Bi{x) = ±\x\t[J_ i (t)-J + i(0], x<0. (50) 

In terms of Airy functions, the solutions (44) and (46) can be written as 

Z+ = '^{C 2 - G\)Ai(x) + ^-(C 2 + Ci)Bi{x) for x > 0, (51) 

Z+ = \{D 2 + D 1 )Ai{x) + ^y(D 2 - Di)Bi(x) for x < 0. (52) 

By matching boundary conditions it is easy to show that the solution corresponding x > and that corresponding 
x < are continuous when C\ = — D\ and C 2 = D 2 . 

As an example of solutions from this region, we choose: a = 0.95, M = 1, m p = 0.17, I = |, m = — ^, and er = 0.21. 
The black hole horizon is at r + = M + V(M 2 - a 2 ) = 1.31225, a c = 0.526316, a s = 0.180987, a 2 = -1.356 and 
A = 0.93 [7]. It is clear that the values of cr c , cr s and a 2 indicate the region is non super-radiant. In Fig. 5a, we show 
the nature of V + and V-, however, while solving, we use the equation containing V + (eq. 24). Unlike the case in the 
previous sub-section, here a 2 is no longer greater than V± at all radii. As a result, k 2 may attain negative values in 
some region. In Fig. 5b, nature of V+ (solid curve), parameter k (dashed curve) and energy E (short-dashed curve) 
are shown. Here, WKB approximation can be applied in regions other than f* <~ —4 to —1 and 1 to 7 where k is close 
to zero and the condition jjf- « k is not satisfied. 

In the region f* = 7 to 1 around the turning point f* = 4.45475 the solutions turns out to be [10], 

Z + = 1.087526Ai(a;) + 0.788968Bi(x). (53) 

Similarly, the solution from f* ~ — 1 to — 4 i.e. around the turning point f* = —2.8053 turns out to be [10], 

Z + = -[1.328096Ai(a;) + 0.774426Bi(x)]. (54) 
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It is to be noted that in the region f* <~ 1 to — 1, even though the potential energy dominates over the total energy, 
WKB approximation method is still valid. Here the solution will take the form ° p J t "* and ^^j— ■ Asymptotic 
values of the instantaneous reflection and the transmission coefficients (which are traditionally known as the 'reflection' 
and 'transmission' coefficients) are obtained from the WKB approximation. This yields the integral constant c as in 
previous case. ^From eq. 37(a-b) local reflection and transmission coefficients are calculated, behaviours of which are 
shown in Fig. 5c. The constants A-h and A + h are calculated as before. Note the decaying nature of the reflection 
coefficient inside the potential barrier. 

Solutions in the Super-Radiant Region 

In this region, the potential diverges at r = \ot\. Here, the barrier height goes up to infinity and then the potential 
changes sign so that its nature changes from repulsive to attractive and vice versa. This is because a < <r s (which 
is the condition for super-radiance) and r + < \a\ [see equations (11) and (13)]. Unlike the previous two cases, the 
relation between r and f* is not single valued. Here, at both r = r + and r = oo, the value of f* = oo. With the 
decrement of r, f* is decreased initially up to r — \a\. Subsequently, starts to rise and at the black hole horizon it 
diverges. Obviously, in this case particles hit the barrier and we can solve the equation following the same methods 
as explained in the previous cases, i.e., eqs. (35) and (35') for the region where the WKB method is valid and eqs. 
(51) and (52) where the WKB method is not valid. 

For illustrative example, here, we choose: a — 0.95, M = 1, m p — 0.105, I = |, m = — ^, and a = 0.105. The 
black hole horizon is located at r+ = M + ^/(M 2 - a 2 ) = 1.31225, and a c = 0.526316, a s = 0.180987, a 2 = -3.62 
and % = 0.97 [7]. Chandrasekhar showed [4] that for integral spin particles this region exhibits super-radiance and 
conjectured that for half-integral spins the super-radiance may be absent. We investigate here if this conjecture is 
valid. 

The behaviour of potentials V+ and V- are shown in Fig. 6a. It is clear that at r — \a\ the potential diverges and 
the nature of the potential is changed from repulsive to attractive (for V-) and vice versa (for V+). Here, we will treat 
the equation with V- as the potential (it is equally easy to do the problem with V + ). We first divide our computations 
into two parts. In the repulsive part of the potential (i.e., when V- > 0), particles come from infinity and most of 
them reflect back from the infinitely high barrier. In the attractive part of the potential (i.e., when V+ < 0), particle 
radiates outwards in the f * co-ordinate (actually, particle goes towards the horizon but due to multivalueness of the 
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radial co-ordinate r* (with respect to r) the horizon is mapped to infinity). 

In Fig. 6b, nature of V-, k and E are shown. The WKB approximation (more precisely IWKB approximation) 
method is valid from infinity to r* ~ 40 since, otherwise, \-§f- « k is not satisfied. In those other regions one has 
to apply a different method (which was also explained in last sub-section) to find solutions. The local values of the 
reflection and transmission coefficients and the wave function of the particle are calculated as in the previous cases. 
Since the matter which tunnels through the infinitely high barrier face infinitely strong attractive field, the possibility 
of extraction of energy would be zero. In Fig. 6c, the variations of local transmission and reflection coefficients are 
shown. The net transmission of the wave through the horizon is non-negative all along and therefore super-radiation 
is absent, although a is less than a s . We believe that the non-existence of super-radiation is due to (r — |a|)~ 3 
variation of the potential near the singular point. Because of the existence of attractive field, the extraction of energy 
is difficult, so the net transmission of the wave through the horizon from oo is always positive. This argument is valid 
for any set of parameters where a < a s . 

IV. CONCLUSION 

In this paper, we studied scattering of massive, spin-half particles from a Kerr black hole, particularly the nature 
of the radial wave functions and the reflection and transmission coefficients. Our main motivation was to give a 
general analytical expression of the solution which can be useful for further study. We showed few illustrative cases 
as examples. We verified that these analytical solutions were indeed correct by explicitly solving the same set of 
equations using quantum mechanical step-potential approach as described in Section III. We classified the entire 
parameter space in terms of the physical and unphysical regions and the physical region was further classified into 
two regions, depending on whether the particle 'hits' the potential barrier or not. Again, the region where particle 
hits on the barrier, is divided into two parts, one is super-radiant region and other is non-super-radiant region. We 
chose one illustrative example in each of the regions. We emphasize that the most 'interesting' region to study would 
be close to m p ~ a. However we pointed out (Fig. 1) that for m p < 0.35, the WKB solutions cannot be trusted, and 
other methods (such as those using Airy functions) must be employed. 

We used the well known WKB approximation method as well as the step-potential method of quantum mechanics 
to obtain the spatial dependence of the coefficients of the wave function. This in turn, allowed us to determine the 
reflection and transmission coefficients and the nature of wave functions. The usual WKB method with constant 
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coefficients and (almost) constant wave number k is successfully applied even when the coefficients and wave number 
are not constant everywhere. Solution from this 'instantaneous' WKB (or IWKB) method agrees fully with that 
obtained from a purely quantum mechanical method where the potential is replaced by a collection of steps. Our 
method of obtaining solutions should be valid for any black hole geometry which are asymptotically flat so that 
radial waves could be used at a large distance. This way we ensure that the analytical solution is close to the exact 
solution. In Region II, in some regions, the WKB method cannot be applied and hence Airy function approach or 
our step-potential approach could be used. 

In the literature, reflection and transmission coefficients are defined at a single point. These definitions are mean- 
ingful only if the potential varies in a small region while studies are made from a large distance from it. In the present 
case, the potential changes over a large distance and we are studying in these regions as well. Although we used the 
words 'reflection' and 'transmission' coefficients, in this paper very loosely, our definitions are very rigorous and well 
defined. These quantities are simply the instantaneous values and in our belief more physical. The problem at hand 
is very similar to the problem of reflection and transmission of acoustic waves from a strucked string of non-constant 
density where reflection and transmission occurs at each point. 

Among other things, we verify Chandrasekhar's conjecture [4] based on the asymptotic solution, that for spin-i 
particle the phenomenon of super-radiance is absent. We believe that this is due to the very way the potential develops 
the singularity at r = \a\. Here I^-(r*) oc (r — H)~ 3 , which results in an attractive potential in some region very 
close to the black hole. In contrast, V-(r*) oc (r — |a|) -4 when electromagnetic and gravitational waves are scattered 
off the black hole [4] does not create an attractive part in the potential and possibly exhibit the phenomenon of 
super-radiance. It is noted that all the cases where potential diverge at r = a (i.e., so called super-radiation cases) 
arise for a < u s with the negative values of azimuthal quantum number (here, m = —1/2) and the positive Kerr 
parameter, a. For positive values of m and positive values of a, potential does not diverge at any point for all values 
of a. If we change the spin orientation of the black hole (negative values of a) and take positive m again divergence of 
the potential will arise. Thus, it seems that the cases with opposite sign of a and m are physically more interesting. 

It is seen that for different physical parameters the solutions are different. The waves scattered off are distinctly 
different in different parameter regions. In a way, therefore, black holes can act as a mass spectrograph! Another 
interesting application of our method would be to study interactions of Hawking radiations in regions just outside the 
horizon. 
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FIGURE CAPTIONS 



Fig. 1: Contours of constant w max — max(p^) are shown to indicate that generally to << 1 and therefore the 
WKB approximation is valid in most of the physical region. Labels indicate values of w max . 

Fig. 2: Behaviour of (a) V + (solid curve) and V- (dashed curve), (b) V + (solid curve), k (dashed curve), total energy E 
(short-dashed curve), (c) local transmission (T, solid curve) and reflection (R, dashed curve) coefficients as functions 
of r*. The parameters are a = 0.5, M = 1, m p = 0.8, I = |, m = — 5 , a = 0.8. 

Fig. 3a: Steps (solid) approximating a potential (dotted), thus reducing the problem to that of a quantum mechanics. 
The parameters are same as in Fig. 2. 

Fig. 3b: Comparison of variation of instantaneous reflection coefficient R and transmission coefficient T with the 
radial coordinate f* using analytical WKB method (solid) and step-potential method (dotted). The parameters are 
same as Fig. 2. 

Fig. 4: Nature of real and imaginary parts of radial wave functions for Case 1. 

Fig. 5: Plots are same as in Fig. 2. The parameters are a — 0.95, M = 1, m p = 0.17, I = \,m = — |, a = 0.21. 
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Fig. 6: Behaviour of (a) V + (solid curve) and V_ (dashed curve), (b) V- (solid curve), k for region where potential 
is positive ( k rep , dashed curve), k for region where potential is negative ( k at t, short-dashed curve), total energy E 
(dotted curve), (c) local transmission (T, solid curve) and reflection (R, dashed curve) coefficients as functions of r*. 
The parameters are a = 0.95, M = 1, m p — 0.105, I = |, m = — |, cr = 0.105. 
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